Optomechanical tailoring of quantum fluctuations 
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We propose the use of feedback mechanism to control the level of quantum noise in a radiation 
field emerging from a pendular Fabry-Perot cavity. It is based on the possibility to perform quantum- 
nondemolition measurements by means of optomechanical coupling. 
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I. INTRODUCTION 



Optomechanical coupling between light and (movable) mirrors gives rise to a variety of intriguing effects. A wide 
class of quantum states coming from such coupling has been investigated in Ref. Jjj. Furthermore, due to the recent 
technological developments in optomechanics, this area is now becoming experimentally accessible Q and quite 
I , interesting. Usually in such a systems, a mirror moves in response to radiation pressure exerted by the light. This 
' leads to an intensity-dependent phase shift for the light field equivalent to an optical Kerr effect [||. Hence, among 
the others, this system shows bistable behaviour Q, and quantum noise reduction effects |]||; moreover, it can be 
used to realize quantum-nondemolition (QND) measurement on the light field Some of these apects have been 
recently studied also in the spatial domain || . 

On the other hand, the use of feedback to control the noise in optical system has been long used. Classical 
theory is often adequate to describe such feedback. However, the advent of nonclassical light (in particular squeezed 
light), required a more accurate treatment. A number of useful approaches exist, using linearized quantum Langevin 
equations [plJlOfl and, more recently, exact quantum trajectories as well [ fllf . In particular, it has been shown that it 
is possible to use QND-mediated feedback to get squeezing, both extracavity and intracavity . 

Along this line, here we investigate the possibility of controlling the light statistics via optomechanical coupling and 
£^ ' a consequent feedback action. 

The paper is organized as follows: in Section II we present the measurement model, in Section III we study the 
light field dynamics, and the spectra are shown in Section IV. Then, in Section V we introduce the feedback action 
and show its influence on the light statistics. Section VI concludes with a numerical result. 



II. THE MEASUREMENT MODEL 

The system we wish to study consists of a Fabry-Perot cavity with a partially transmitting mirror on one side and 
a perfectly reflecting, and oscillating, mirror on the other side. We imagine the latter coated on the surface of a 
piezoelectric crystal and then we shall adopt the model of Ref. |l3| . In such a case the variation of the crystal length 
can be measured by an electric circuit. Hence, in reality, the system is composed of three different coupled subsystems, 
the light, the crystal and the electric circuit. The light field is coupled to the crystal length by the radiation pressure 
force, while the crystal and the electric circuit are coupled by piezoelectric effect. 

Concerning the light field, the system is equivalent to the usual Fabry-Perot cavity except that the cavity detuning 
may vary because of the crystal elasticity. The dynamics of the intracavity field a in a rotating frame is described by 

^ = -[ 7 -i T -i*(t)]a(t) + ^a in (t), (1) 

where r is the cavity round trip time and 7 is the amplitude damping rate of the cavity (equal to the transmissivity of 
the fixed mirror divided by r). Furthermore, ai n is the incoming field, and it is related to the intracavity and outgoing 
ones by the input-output relation 



flout (*) = y/2-ya(t) - a ia (t) . (2) 
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The cavity detuning \&/t depends linearly on the variation x(t) of the crystal length 



*(f) = tA + 2k x(t) , (3) 

and couples the mechanical motion of the crystal to the field. Here, fco = ojq/c is the incoming field wave number, 
while A is the cavity detuning in the absence of the mirror's motion. 

The mechanical response of the crystal is described by the variation x(t) of the crystal length. For small displace- 
ments, the linear response theory can be invoked 

x(u) = xoM [Fr(cj) + CQM + Ft(u)] , (4) 
where Xo(w) 1S the mechanical susceptibility 

1/m 

XoM = — = : , (5) 

uj^ - lo z - iuj m 

of the (movable mirror) crystal which has mass m, oscillates at frequency uj m , and it is damped at rate j m . 
The first force acting on the mirror, Fr, is the radiation pressure force given by 

F R = 2hk a)a/T . (6) 

This has the simple interpretation of the momentum change of a photon under reflection, times the number of photons, 
divided by the round-trip time. 

The electric circuit induces the presence of an electric charge Q on one side of the crystal and -Q on the other side 
(the current / being the time derivative of Q, i.e. I(ui) = —iu)Q{uj)). These charges generate a variation of the crystal 
length proportional to Q. This effect is equivalent to a force applied on both sides of the crystal, with opposite sign 



and absolute value (Q, where ( is a constant depending on the crystal characteristics |15|. 

Finally, Ft corresponds to the Langevin force describing the coupling of the mechanical motion with a thermal 
bath. Its spectrum is given by 

S Ft = 2mj m k B T, (7) 

where T is the temperature of the thermal bath and ks the Boltzmann constant. This relation is valid when T is such 
that ksT ^> hui m . Otherwise, one should bring some corrections to the dynamics of the mechanical oscillator ]l7| . 

Now, the piezoelectric crystal is equivalent to a capacitance C and a voltage generator in series. A length variation 
induces a polarization inside the crystal which generates a potential difference between the two sides of the crystal. 
This voltage is proportional to the length variation (V — —(x) and thus, a measurement of this voltage allows one 
to monitor the crystal motion, and to measure the amplitude quadrature of the intracavity field (or equivalcntly the 
radiation pressure force) . 

The measured signal is the current I out going out of the system by a coaxial line; a counterpropagating current P n 
enters the system by the same line. This latter corresponds to the Nyquist noise of the characteristic resistance R of 
the line. The input-output relation for the current is jD| 

i = i in + r ut , (8) 

and the noise spectrum of the current I m entering the system is given by 

S„- k -§. (9) 

The voltage at the entrance of the circuit is related to the currents by two equations, one for the line and one for the 
circuit 

V(u)=R[l i »(u))-I™ t (Lj)} , (10) 
V{uj) = Z (w)I(w) - (x(u) , (11) 

where Zq is the (purely imaginary) impedance of the circuit including the capacitance C of the crystal, and —C,x is 
the voltage generated by the piezoelectric effect. 

From Eqs.(||), (poh , (O), and with the aid of Eq.fl) it is possible to get the output, measured current as 
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t = r-z(u) . CxoM [F ( } F ( }] (12) 

V ; R+Z{uj) V 7 i? + Z(w) L -RV > V ) 



where Z(cu) is an effective impedance 



c 2 

Z(w) = Z (u) - i— Xo(w) • (13) 



We see that the measured current provides information about the radiation pressure force. 



III. THE LIGHT FIELD DYNAMICS 

Following Ref. p3| , to study the light field dynamics we exploit the linearisation procedure. The equation for the 
mean field, which follows from Eq. (uy, is 



= - [7 - ty/r] a + y/2j a in , (14) 
where ip represents the steady state detuning. The solution reads 

a = r-7-ain ■ (15) 

7 — I1p/T 

For the sake of simplicity, we shall assume, from now on, a- m to be real. After linearisation, we get the equation for 
the quantum fluctuations 

j 5a(t) = - [7- ir _ V] Sa(t) + iaS^(t)/r + ^ Sa in (t) . (16) 
Now, by defining the generic quadrature 

X e = (ae- w + a)e ie ) , (17) 

the radiation pressure force (|J) becomes 

F R = 2hk Q T- 1 [\a\ 2 +a5a) +a*5a] = 2hk T~ 1 \a\ (\a\+5X avga ) , (18) 

so, it is proportional to the fluctuations of a quadrature determined by the stationary value of the intracavity field. 
The equation of motion (|l6|), can be rewritten as 

— 6X 9 (t) = -7 ( 5A e (t)-r- 1 V"5A e+V2 (t) + 2|a|r- 1 sin(0-arga) 5*(t) + ^/¥(8Xf{t) , (19) 

where the spectrum of the input fluctuations is <Sx in = 1 f° r coherent light. We immediately see that the measured 
quadrature 9 — arga is not affected by the measurement process (i.e. being a QND observable) provided to have 
ip = 0. So, from now on we assume such a condition, which also implies a to be real (i.e. arga = 0). Thus, in the 
frequency domain, the solution of Eq. ( jigj ) reads 

5X 9 {u) = 2asin * 8*{w) + T^h ■ (20) 

(7 — IU1)T (7 — IU>) 

Any real measurement device must give a read-out and in this case it is the current I out . Hence the fluctuations of 
this current is directly related to the fluctuations of the quadrature Xq. Specifically, by using Eqs. (p"2]), and (p8|), we 
have 

5I ou \lu) = a ia g(u>)8X™(cj) + I T (u>) , (21) 
where the gain Q, and the thermal current 7y, are given by 

g(u) = 2Kk ( 1 (22) 

(7 — VjJ)T K + Z(uj) 
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R + Z(oj) R+ Z(uj) 

In the case of equal mechanical and electric temperatures, the spectrum of the thermal current can be easily calculated 
by means of Eqs.(||), (Q), and (|l3|), and the result 

as required for thermal equilibrium. 



IV. THE LIGHT SPECTRUM 



In order to derive the light spectra we first calculate the spectrum and the correlations x ^ , C* i T |19| . To 
this end, we remind that the current / is the time derivative of the charge Q (i.e. I(w) = —iujQ{lo)) appearing in 
Eq.(^), and then we eliminate the electric variable in Eqs.([|), and (|IT|), obtaining 

.C 2i? 



= xH 



F r {lo) + F t (lu) + i 



lu R + Zq(u>) 



where 



1 



1 



• C 2 



f 



(25) 



(26) 



2R 



xM XoM u> R + Z (u)) 

This means that the mechanical response to the radiation pressure is modified by the coupling. As a consequence it 
will be 

6^(u>) = 2fc x(w 
Now by considering Eqs.(|27]), (§|), 



2hk T~ 1 \a\5X {uj) + F t {uj) + i- 



uj R + Zo(uj) 

(p0|), (||), and (f7j) we can calculate the desired quantities 



Ah 2 kl\a\ 2 
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c 2 



2R 



T~2~i — 2T~2 + 2m7 m fc B T . . .„ 



k B T 



-9 XT 



Ahk 2 a X (co) 



(7 — Iuj)t 



and 



C*/ T = 2/c x(w) 



Cx5H 



.c 



I R-Z*(u) 



_R + Z*(lu) ' Wj R+Zo(w)i? + Z*(w) 
Now, from Eq. ( po| ) the intracavity spectrum can be written as 



k B T. 



(27) 

(28) 
(29) 

(30) 



S 



x 



Aa 2 sin 2 9 
(7 2 + cj 2 )t 2 



<S«p(a;) + 
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7 2 + a; 2 



2 v / 27asin( 



2Re 



( 7 2 +w 2 )r 

Finally, the output spectrum can be derived by using the input-output relation . The result is 

87a 2 sin 2 9 



(31) 



(7 2 + oj 2 )t 2 
2y/2ja sin 9 



(7 2 + uj 2 )t 



2Re{( 7 + iu;)C* X i»} 



(32) 



It is immediately apparent that the output spectrum (32) remains equal to one for 9 = [a fact obvious from Eq. ( J2C 
and used in Eq. (28)], while it takes its maximum value for 9 = n/2. This is due to the back action of the mirror noise, 
which in the regime we consider is dominated by thermal noise. Furthermore, calculating Eq.(^2|) we neglected the 
thermal photons generated by the mirror motion. This is reasonable at optical frequencies (ojq) since the mechanical 
frequency (cu m ) is much smaller ||. 

Notice, that in addition to Ref. jl3| we explicitely calculated the spectrum for the detuning and for a generic 
quadrature. 
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V. THE FEEDBACK LOOP 



We now assume that the cavity mode is controlled by modulating the input field linearly with the measured current. 
This can be achieved by using electro-optic modulators [20| . To model this feedback mechanism, we use the quantum 
Langevin approach rather than the quantum trajectory approach [pl| . Moreover, to simplify the problem, we assume 
that the feedback is Markovian, that is the response function of the feedback loop is approximately flat from zero to 
a frequency much larger than that of the mirror or of the electric circuit. This is not unreasonable, as electronic and 
electro-optic devices can respond on time scales much shorter than 10 _6 s. Hence, the feedback can be modeled by 
changing the input beam according to the following rule 



a in (t) =a in + Sa in (t) a in (t) = a in + 6a- m (t) - XSI out {t) . 

In terms of the quadrature fluctuation operator 

5Xf(t) -» 5Xf(t) = 5Xf(t) - (\e- te + AV e ) SI out (t) . 



(33) 



(34) 



Here, the value of A is related to the practical way of realize the feedback action, i.e., to the electro-optics modulators 
capabilities. 

With the choice of A real, the quadrature Eq.(|20|) is modified as follows 



SX e (u) 



2a sin 6 

(7 — ilo)t 
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-5Xf{u) 



2A ^l[cos ( 

7 — ILU 



(7 - iu) 

/(a;) sin0] [a in g(uj)SX^(cu) + I t (lu)] , 



where 



f(u) = 8hk 2 



2 a 2 x(w) 



(7 — zw)r 2 



(35) 



(36) 



is the feedback term coming from the modified detuning fluctuations, i.e. the term proportional to A when the 
replacement (^4|) is made in Eq.(p7[). 
The intracavity spectrum then becomes 



4A 

-4A 
+ 4A 2 



20 



7 2 + id 2 



Re {a- m g*{uj)e~ iB [cos 6 + sin (9]} 

1 ' " Re { [sin (26) + 2f* (u) sin 2 6] \a in G* (wJC* x «» (u) + It (w)1 } 



( 7 2 + uj 2 )t 
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7 2 + LO 2 



«fn I G(L0) I 2 + 



2R 



|cos0 + /(w)sin0|' 



In view of Eqs.(34), 

SX2 ut = SX2 ut - 2\ 



, the input-output relation (g) gives 

'j + iu 27/(0;) . . 
cos6» + u \ ' sing 

7 — iuj 7 — ILU 

Finally, the output spectrum in the presence of feedback action is 

cos 6 H 



[a in g(u)5X™(uj) + I T (u)] . 



(37) 



(38) 



Sxovt = 5 x; »t - 4ARe <^ a- m Q*(u>)e 



■ sin ( 



4ARc 
4A 2 



27a 



ol I G{u>) 



7 — ilo 

sin (26) 47/*(w)sin 2fl1 
7 + iu> 

12 , k sT 



7 2 + LU 2 



2R 



1 + a 2 7 /(w) . 
cos 6 H sin 



a-mG* (oS) C w X i„ (uj) + Cq, It (u>) 

2 



7 — ILU 



7 — no 



(39) 



For the quadrature X the above expression reduces to 
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s 



y,,„t 



1 - 2a in XG(u J ) 
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(40) 



For any given frequency u), it is possible to choose some value of feedback strength A such that this expression 
will be less than one, providing the real part of Q(oj) is positive. That is, the feedback enables one to produce a 
sub-shot-noise output. The optimum value for the feedback parameter results 



A pt — 



QinRc [G(uj)} 



(41) 



For the amount of squeezing to be significant it is necessary for the measurement gain to be mostly positive and also 
large compared to the thermal noise: 



a in \G(uj)\ > 



2R 



(42) 



The same conditions hold for the intracavity noise as can be seen from Eq.(E 



VI. CONCLUSION 



The main result, obtained in the preceding section, is that using an electro-optic feedback mechanism it is possible 
to reduce the quantum fluctuations of a field quadrature below the shot-noise limit. The amount of noise reduction 
depends on the value of the feedback parameter A as well as on the frequency oj determining the gain in the measure- 
ment process. The best squeezing is obtained in the amplitude quadrature (with respect to the input field) for which 
a perfect QND measurement is achievable. On the other hand a rotation of the ellipse of fluctuations [£2| could be 
devised at the output. It is also worth noting that in the model we have elaborated, the added thermal noise (Eq.(|40|)) 
can be controlled by means of the electric components (e.g. the resistance), contrary to the case of Refs. where 
thermal fluctuations lead to anavoidable detrimental effects. 

To better show the potentialities of our model, we consider a realistic experimental situation. First of all, to enhance 
the measurement efficiency, we imagine we have a resonant electric circuit. With self-induction of impedance L the 
electric resonance frequency oj e is given by 



LCtoi = 1 , 

and we choose this to be almost equal to uj m . The noncoupled electric impedance Zq is 



(43) 



Z (u) 



1 



Luj 



(44) 



Then, it is possible to define the electric damping and the piezoelectric coupling frequency 



7e = R/L , n e = evej 



m . 



Now, we take for granted the values of parameters in Ref. (l3|. The mirror parameters are: 
value should be intended as an effective value coming from the acoustic modes of the mirror [E 



= 10 b s" 



quality factor Q m = Um/jm = 10 
-l 



and piezoelectric coupling frequency f2 e 

"He 



(45) 

mass m = 10~ 3 g (this 
| ) , resonance frequency 
« 10 3 s _1 . The cavity 



parameters are: bandwidth 7 = 10 6 s _1 , round trip time r = 10 _il s, field wavelength 0.5 /im, and incident power 
P m = hu}oaf n = 100 mW. The transmissivity of the fixed mirror should be the smallest possible since the gain factor 
( p2| ) is inversely proportional to it. Finally the electric quality factor Q e = ui e /-f e w 10 6 . 

In Fig.l we display the spectrum of Eq.(40) for different values of temperature. We see that the noise reduction 
occours in proximity of the mechanical resonace frequency, where the response of the meter is maximum. It is worth 
noting that the squeezing bandwidth is approximately proportional to the meter bandwidth j m , and whenever it is 
increased the squeezing phenomenon becomes more sensible to the thermal noise since the mechanical quality factor 
Q m decreases. 

Finally, we would briefly discuss about the quantum mechanical consistency of our treatment. For the QND 
variable Xq, and its conjugate the commutation relations are preserved from the input to the output. The same 
is not true in general for Xg and Xg +7r /2 with 0^=0, due to the fact that the mirror introduces beside a frequency 



G 



dependent phase shift, an additional frequency-dependent damping Hence, the canonical commutation relation 

[Xg nt (t), X°%(t')] = 2iS(t - £') holds, and it leads to the uncertainty relation 

5 X out (w)5 X o„t (w) > 1 , (46) 

for all oj. In reality, since there is a second output from our system, the current 7 out , which measures Xy ut , we can 
demand a stronger inequality, 

S x ™t\i<m% H<S X out (w) > 1 . (47) 

Here the first quantity is the spectrum of fluctuations in the 9 = quadrature conditioned upon the measured current 
I out . For Gaussian statistics, which hold in the linearized regime we are considering, this conditioned spectrum can 
be evaluated as 

5 X out| JO ut(a;) = S X out(u;) - 

We have checked the relation (|47| ) numerically for the regimes of interest, both with and without feedback, and find 
that it is always satisfied. 

In conclusion, we have presented a realistic model, based on optomechanical coupling, for an active control of 
the quantum noise of a radiation field. The theory has been developed by using a semiclassical approach, which 
nevertheless guarantes quantum mechanical consistency for realistic situation. Of course in a real situation minor 
detrimental effects could appear; we mention for example the non perfect reflectivity of the rear mirror (which decreases 
the gain factor), or the absorption and diffraction losses in the input coupler (which introduce excess noise). However, 
we think that the up-to-date technology is mature enough to envisage experimentally an optomechanical tailoring of 
quantum fluctuations; this could be useful in cavity based experiments, or to produce nonclassical light; among other 
possible interesting applications we cite the reduction of the radiation pressure noise in gravitational interferometers 




ut (Ll>) I Slant (Ll>) 
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